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CRITICAL SURFACE OF THE 
HEXAGONAL POLYGON MODEL 


GEOFFREY R. GRIMMETT AND ZHONGYANG LI 

Abstract. The hexagonal polygon model arises in a natural way via a transfor¬ 
mation of the 1-2 model on the hexagonal lattice, and it is related to the high 
temperature expansion of the Ising model. There are three types of edge, and 
three corresponding parameters a, /3,7 > 0. By studying the long-range order of a 
certain two-edge correlation function, it is shown that the parameter space (0, 00)^ 
may be divided into subcritical and supercritical regions, separated by critical sur¬ 
faces satisfying an explicitly known formula. This result complements earlier work 
on the Ising model and the 1-2 model. The proof uses the Pfaflian representation 
of Fisher, Kasteleyn, and Temperley for the counts of dimers on planar graphs. 


1. Introduction 

The polygon model studied here is a process of statistical mechanics on the space 
of disjoint unions of closed loops on hnite subsets of the hexagonal lattice El with 
toroidal boundary conditions. It arises naturally in the study of the 1-2 model, and 
indeed the main result of the current paper is complementary to the exact calculation 
of the critical surface of the 1-2 model reported in [5, 6] (to which the reader is referred 
for background and current theory of the 1-2 model). The polygon model may in 
addition be viewed as an asymmetric version of the 0{n) model with n = 1 (see [2] 
for a recent reference to the 0{n) model). 

Let G = (V, E) be a finite subgraph of El. The configuration space Eg of the 
polygon model is the set of all subsets S of E such that every vertex in V is incident 
to an even number of members of S. The measure of the model is a three-parameter 
product probability measure conditioned on belonging to Eg, in which the three 
parameters are associated with the three classes of edge (see Figure 2.1). 

This model may be regarded as the high temperature expansion of a certain in- 
homogenous Ising model on the hexagonal lattice. The latter is a special case of the 
general eight-vertex model of Lin and Wu [16] (see also [26]). Whereas Lin and Wu 
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concentrated on a mapping between their eight-vertex model and a generalized Ising 
model, the current paper utilizes the additional symmetries of the polymer model 
to calculate in closed form the equation of the critical surface for a given choice of 
order parameter. The parameter space of the polymer model extends beyond the set 
of parameter values corresponding to the classical Ising model, and thus our overall 
results do not appear to follow from classical facts (see also Remarks 2.2 and 2.4). 

The order parameter used in this paper is the one that corresponds to the two- 
point correlation function of the Ising model, namely, the ratio Zc^e^f jZ q-, where 
ZG,e*^f is the partition function for conhgurations that include a path between two 
edges e, /, and Zq is the usual partition function. 

Here is an overview of the methods used in this paper. The polymer model may 
be transformed into a dimer model on an associated graph, and the above ratio may 
be expressed in terms of the ratio of certain counts of dimer conhgurations. The last 
may be expressed (by classical results of Kasteleyn [7, 8], Fisher [3], and Temper ley 
and Fisher [20]) as Pfafhans of certain antisymmetric matrices. The squares of these 
Pfafhans are determinants, and these converge as G ^ I® fo the determinants of 
inhnite block Toeplitz matrices. Using results of Widom [23, 24] and others, these 
limits are analytic functions of the parameters except for certain parameter values 
determined by the spectral curve of the dimer model. This spectral curve has an 
explicit representation, and this enables a computation of the critical surface of the 
polygon model upon which the limiting order parameter is non-analytic. 

More specihcally, the parameter space (0, cx))^ may be partitioned into two regions, 
termed the supercritical and subcritical phases. The order parameter displays long- 
range order in the supercritical phase, but not in the subcritical phase. 

The results of the current paper bear resemblance to earlier results of [6], in which 
the same authors determine the critical surface of the 1-2 model. The outline shape 
of the main proof (of Theorem 2.5) is similar to that of the corresponding result of 
[6]. In contrast, neither result seems to imply the other, and the dimer correspon¬ 
dence and associated calculations of the current paper are based on a different dimer 
representation from that of [6]. 

The characteristics of the hexagonal lattice that are special for this work include 
the properties of trivalence, planarity, and support of a action. It may be possible 
to extend the results to other such graphs, such as the Archimedean (3,12^) lattice, 
and the square/octagon (4,8^) lattice. 

This article is organized as follows. The polygon model is dehned in Section 2, 
and the main Theorem 2.5 is given in Section 2.3. The relationship between the 
polygon model and the 1-2 model, the Ising model, and the dimer model is explained 
in Section 3. The characteristic polynomial of the corresponding dimer model is 
calculated in Section 3.5, and Theorem 2.5 is proved in Section 4. 
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2. The polygon model 

We begin with a description of the polygon model. Its relationship to the 1-2 
model is explained in Section 3.1. The main result (Theorem 2.5) is given in Section 

2.3. 

2.1. Definition of the polygon model. Let the graph G = {V,E) be a hnite 
connected subgraph of the hexagonal lattice El = (V, E), suitably embedded in 
as in Figure 2.1. The embedding of El is chosen in such a way that each edge may 
be described by one of the three directions: horizontal, NW, or NE. (Later we shall 
consider a hnite box with toroidal boundary conditions.) Horizontal edges are said 
to be of type a, and NW edges (respectively, NE edges) type b (respectively, type c), 
as illustrated in Figure 2.1. Note that El is a bipartite graph, and we call the two 
classes of vertices black and white. 

Let n be the product space H = {0,1}^. The sample space of the polygon model 
is the subset nP°h = nP°h(G) C H containing all tt = (tTe : e G E) G H such that 

(2-1) E TTe is either 0 or 2, v G V. 

eBv 

Each 71 G nP°h j^iay be considered as a union of vertex-disjoint cycles of G, together 
with isolated vertices. We identify tt G H with the set {e G E : TTe = 1} of ‘open’ 
edges under tt. Thus (2.1) requires that every vertex is incident to an even number 
of open edges. 



Figure 2.1. An embedding of the hexagonal lattice, with the edge- 
types marked. 

Let ea, Cb, Cc 7 ^ 0. To the conhguration vr G -yyg assign the weight 

(2.2) u;(7r) = 

where 7r(s) is the set of open s-type edges of vr. The weight function w gives rise to 
the partition function 

(2.3) Zo{P)= E 

7renp°>y 
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This, in turn, gives rise to a probability measure on nP°h given by 

(2.4) PG(,r) = 

The measure may be viewed as a product measure conditioned on the outcome 
lying in 11 ^°^ yyg concentrate here on an order parameter to be given next. 

It is convenient to view the polygon model as a model on half-edges. To this end, 
let AG = {AV^ AE) be the graph derived from G = (V, E) by adding a vertex at the 
midpoint of each edge in E. Let ME = {Me : e G E} be the set of such midpoints, 
and AV = V U MV. The edges AE are precisely the half-edges of E, each being of 
the form {v, Me) for some v ^ V and incident edge e & E. A polygon conhguration 
on G induces a polygon conhguration on AG, which may be described as a subset of 
AE with the property that every vertex in AV has even degree. For an a-type edge 
e E E, the two half-edges of e are assigned weight (and similarly for b- and c-type 
edges). The weight function w of (2.2) may now be expressed as 

(2.5) M;(7r) = tteW^^^AG). 

We introduce next the order parameter of the polygon model. Let e, f E ME he 
distinct midpoints of AG, and let Lie,/ be the subset of all vr G {0, 1 }"^® such that: 

(i) every v E AV with v ^ e, f is incident to an even number of open half-edges, and 

(ii) the midpoints of e and / are incident to exactly one open half-edge. We dehne 
the order parameter as 

( 2 . 6 ) = 

where 

(2.7) := 

TTSlIe,/ 

Remark 2.1 (Notation). We write for the parameter of an s-type edge, and eg for 
that of edge g. For conciseness of notation, we shall later work with the parameters 

2 Q 2 2 

Ot . 6^, p . 6^, p . €g, 

and the main result. Theorem 2.5, is expressed in terms of these new variables. The 
‘sguared’ variables ef are introduced to permit use of the ‘unsguared’ signed variables 
Cs in the definition of the polymer model on AG. 

The weight functions of (2.2) and (2.5) are unchanged under the sign change 
Cs —t —Cs for s E {a,b,c}. Similarly, if the edges e and f have the same type, 
then, for vr G hl^j, the weight w{n) of (2.7) is unchanged under such sign changes. 
Therefore, if e and f have the same type, the order parameter Mde, f) is independent 
of the sign of the Cg, and may be considered as a function of the parameters a, (3, 7 . 




CRITICAL SURFACE OF THE HEXAGONAL POLYGON MODEL 


5 


Remark 2.2 (High temperature expansion). If G (0,1)^, the polygon model 

with weight function (2.5) is immediately recognized as the high temperature expan¬ 
sion of an inhomogeneous Ising model on AG in which the edge-interaction Jg of an 
s-type half-edge satisfies tanhJ^ = le^]. Under this condition, the order parameter 
f) of (2.6) is simply a two-point correlation function of the Ising model (see 
Lemma 3.2). If the |es| are sufficiently small, this Ising model is a high-temperature 
model, whence Mc^e, f) tends to zero in the double limit as G f EI„ and \e — f\ —)■ oo, 
in that order. It may in fact be shown (by results of [12, 14, 17] or otherwise) that 
this Ising model has critical surface given by the eguation aft + /97 + 7 a = 1. 

See [1, p. 75] and [18, 22 ] for accounts of the high temperature expansion, and [4] 
for a recent related paper. The above Ising model may be viewed as a special case of 
the eight-vertex model of Lin and Wu [16]. It is studied further in [ 6 , Sect. 4]. 

2.2. The toroidal hexagonal lattice. We will work mostly with a hnite subgraph 
of El subject to toroidal boundary conditions. Let n > 1, and let ti, T 2 be the two 
shifts of El, illustrated in Figure 2.2, that map an elementary hexagon to the next 
hexagon in the given directions. The pair (ti, T 2 ) generates a 1? action on El, and we 
write EI„ = (l/„, En) for the quotient graph of El under the subgroup of generated 
by the powers rf and r^. The resulting EI„ is illustrated in Figure 2.2, and may be 
viewed as a hnite subgraph of El subject to toroidal boundary conditions. We write 
M„(e,/) :=MH„(e,/). 



Figure 2.2. The graph EI„ is an n x n ‘diamond’ wrapped onto a 
torus, as illustrated here with n = 4. 

As in Remark 2.1, let 

(2.8) a = el, /3 = el 7 = el 
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Theorem 2.3. Let e, f be distinct edges o/EI„. The order parameter Mn{e, f) = 
/) is invariant under the change of variables (a,/ 3 , 7 ) h-)■ (a, 7 “^), 
and similarly under the other two changes of variables in which exactly two of the 
parameters a, j3, 7 are replaced by their reciprocals. 

Proof. Let tt G be a polygon configuration on AEI„, and let vr' be obtained from 
TT by 

] 7r(e) if e has type a, 

= otherwise. 

Since tt' is obtained from tt by adding, modulo 2 , a collection of edges that induce 
an even subgraph of ]HI„, we have that vr' G pet be the weight of tt 

as in (2.5), with a, (3, 7 given by (2.8). Then 

(2.10) w“’^’^(7r) = (/?#V")w^“’^“'’^"'(7r'), 

where ffs is the number of s-type edges in EI^. Similarly, if e 7 ^ /, then tt' G Ilej 
and ( 2 . 10 ) holds. 

By (2.9)-(2.10), is unchanged under the map (a,/?, 7 ) i-A (a, 7 “^). The 
same proof is valid for the other two cases. □ 

Remark 2.4. Recall Remark 2.2, where it is noted that the polymer model is the high 
temperature expansion of a solvable Ising model when (a,/?, 7 ) G (0,1)^. By Theorem 
2.3, this results in a fairly complete picture of the behaviour o/hm„^oo Mn{e, f) when 
either none or exactly two of the three parameters lie in (l,oo). In contrast, the 
dimer-based methods of the current work permit an analysis for all triples {a, 13, 7 ) G 
(0, 00 )^, but at the price of assuming that e, f satisfy the conditions of the forthcoming 
Theorem 2.5. 


2.3. Main result. Let e = {x,y) denote the edge e G E of the hexagonal lattice 
El = (V, E) with endpoints x, y. We shall make use of a measure of distance \e — f\ 
between edges e and / which, for definiteness, we take to be the Euclidean distance 
between the midpoints of e and /, with El embedded in in the manner of Figure 
2.2 with unit edge-lengths. Our main theorem is as follows. 


Theorem 2.5. Let e, / G E be NW edges such that: 

there exists a path I = f) of AM.n from Me to M f, 
using only horizontal and NW half-edges. 

Let e<j 7 ^ 0 for s = a,b, c, so that a, /d, 7 > 0, and let 


( 2 . 11 ) 


^ 

0 -|- P 

where 72 is interpreted as 00 if a = (3. 


72 = 


1 -|- oi(3 


a 
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(a) The limit M(e, f)^ = limn^.oo Mn{e, /)^ exists for 7 7 ^ 71 , 72 . 

(b) Supercritical case. Let i?sup be the set of all (a,/?, 7) G (0, 00)^ satisfying 


1 — a(3 
a ft 


< 7 < 


1 ~t“ cx.ft 
a — ft 


The limit A(a, /?, 7) := lim|e_/|^.oo Af(e, /)^ exists on i?sup, Qiic? satisfies A > 0 
except possibly on some nowhere dense subset. 

(c) Subcritical case. Let Rguh be the set of all (a, ft, 7) G ( 0 , 00)^ satisfying 


either 7 < 


1 — aft 
a ft 


or 7 > 


1 + aft 
a — ft 


The limit A{a,ft,j) exists on i?sub (ind satisfies A = 0. 


The function A has a singularity when crossing between the subcritical and super¬ 
critical regions. A brief explanation of the regions i?sub and Rsup follows. It turns 
out that the process is ‘critical’ if and only if the spectral curve (see Section 3.5) 
of the corresponding dimer model intersects the unit torus. This occurs if and only 
the parameter-vector (a, ft, 7 ) is a root of the equation aft -f /37 -|- 7 a = 1 or of any 
one of the three equations obtained from this equation by the changes of variables 
of Theorem 2.3. See Proposition 3.4. 

Assumption (2.11), as illustrated in Figure 2.3, is key to the method of proof, and 
we present no results in the absence of this condition. Thus, Theorem 2.5 is not 
of itself a complete picture of the location of critical phenomena. For the related 
1-2 model, certain further information about the limits corresponding to Theorem 
2.5(b,c) may be derived as described in [ 6 ] and [17], and we do not explore that 
here, beyond saying that it includes information on the rates of convergence, and on 
correlations unconstrained by condition ( 2 . 11 ). 



Figure 2.3. A path £ of NW and horizontal edges connecting the 
midpoints of e and /. 

By Remark 2.1, it will suffice to prove Theorem 2.5 subject to the assumption that 
65 > 0 for s = a, h, c. 






















GEOFFREY R. GRIMMETT AND ZHONGYANG LI 


3. The 1-2 and dimer models 

We summarize next the relations between the polygon and the 1-2 and dimer 
models. 

3.1. The 1-2 model. A 1-2 configuration on the toroidal graph EI^ = {Vn,En) is a 
subset F En such that every n G W is incident to either one or two members of 
F. The subset F may be expressed as a vector in the space = { — 1,1}^" where 
— 1 represents an absent edge and 1 a present edge. (It will be convenient later to 
use the space rather than the more natural n„ = {0,1}'®'".) Thus the space of 
1-2 conhgurations may be viewed as the subset of containing all vectors a such 
that 

^7reG{l,2}, veVn, 

eBv 

where Tie = |(1 + <^e)- 

The hexagonal lattice El is bipartite, and we colour the two vertex-classes black 
and white. Let a, 6 , c > 0 be such that {a,b,c) 7 ^ (0,0,0), and associate these three 
parameters with the edges as in Figure 2.1. For cr G and v G W, let crl^, be the 
sub-conhguration of a on the three edges incident to v, and assign weights tc(cT|^) to 
the cr^ as in Figure 3.1. We observe the states cTe„ Ce^^, cie^^^, where is the edge 
of type s incident to v. The corresponding signature is the word Ti{ey^c)T^{e^,b)T^{^v,a) 
of length 3. The signature of v is given as in Figure 3.1, together with the local 
weight ry(cT|t,) associated with each of the six possible signatures. 



001 , a no, a 010, b 101, b 



100, c oil, c 



001 , a 110, a 



010 , b 


/ 

\ \ 

\ \ 

\ \ 

• • 

101, b 100, c on, c 




Figure 3 . 1 . The six possible local conhgurations a\y at a vertex v 
in the two cases of black and white vertices of EI„ (see the upper and 
lower hgures, respectively). The signature of each is given, and also 
the local weight w{(T\f) associated with each local conhguration. 








Let 

(3.1) 
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w(o-) = 0-e 

vev 

and 

(3.2) Zn= w{a). 

This gives rise to the probability measure 

(3.3) /i„((T) = -^w{a), a e S„. 

We write {X)n for the expectation of the random variable X with respect to 

The 1-2 model was introduced by Schwartz and Bruck [19] in a calculation of the 
capacity of a certain constrained coding system. It has been studied by Li [13, 15], 
and more recently by Grimmett and Li [6]. See [5] for a review. 


3.2. The 1-2 model as a polygon model. By [6, Prop. 4.1], the partition function 
Zn of the 1-2 model with parameters a, b, c on IHI,^ satishes 

Zn = (|(a -I- 5 -I- c))' "^Z'^, 

where 


(3.4) z; = E n y ^l)(J y BCTydCTy^^Q H" C CT V V 

(tGS uGVn 

denotes the state of the s-type edge incident to n G G, and 

a—b—c b—a—c c—a—b 

-, ID = -, G = -. 

o-|-6-|-c Qj -\- b -\- c fl-|-6-|-c 

Each e = {u,v) G En contributes twice to the product in (3.4), in the forms au,s 
and for some s G {a, 5, c}. We write ae for this common value, and we expand 

(3.4) to obtain a polynomial in the variables cXe. In summing over a G S„, a term 
disappears if it contains some with odd degree. Therefore, in each monomial 
M{a) of the resulting polynomial, every ae has even degree, that is, degree either 
0 or 2. With the monomial M we associate the set ttm of edges e for which the 
degree of Ue is 2. By examination of (3.4) or otherwise, we may see that ttm is 
a polygon conhguration in EI^, which is to say that the graph (14 , tm) comprises 
vertex-disjoint circuits (that is, closed paths that revisit no vertex) and isolated 
vertices. Indeed, there is a one-to-one correspondence between monomials M and 
polygon conhgurations tt. The corresponding polygon partition function is given at 
(2.3) where the weights e^, e?,, Cc satisfy 

(3,6) 





B, = G, 
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which is to say that 
(3.7) 


. BC ^ AC . AB 

a - C- 


Note that these squares may be negative, whence the corresponding e^, e^, ec are 
either real or purely imaginary. 

The relationship between and the parameters a, b, c is given in the following 
elementary lemma, the proof of which is omitted. 


Lemma 3.1. Let a > b > c > 0, and let be given by (3.5)-(3.7). 

(a) Let a < b + c. Then Ca, e?,, Cc are purely imaginary, and moreover 

(i) if a? < 6 ^ + C, then 0 < \ea\ < 1 , 0 < \eh\ < 1 , 0 < led < 1 , 

(ii) if = b‘^ + (?, then |ea| = 1 , 0 < |efe| < 1 , 0 < |ed < 1, 

(hi) if of > 6 ^ + C, then jcaj > 1 , 0 < jefej < 1 , 0 < jecj < 1 . 

(b) If a = b + c, then |ea| = oo, = ec = 0. 

(c) //a > 6 + c, then ea,eb,ec are real, and moreover |ea| > 1 , 0 < |efc| < 1 , 

0 < |ed < 1. 

Equations (3.6)~(3.7) express the eg in terms of A, B, C. Conversely, for given 
real e^ 7 ^ 0, it will be useful later to dehne A, B, C by (3.6), even when there is no 
corresponding 1-2 model. 


3.3. Two-edge correlation in the 1-2 model. Consider the 1-2 model on ]HI„ with 
parameters a, b, c, and specihcally the two-edge correlation {ae(Jf)n where e, f & En 
are distinct. 

We multiply through (3.4) by CeCTf and expand in monomials. This amounts to 
expanding (3.4) and retaining those monomials M in which every ag has even degree 
except (Je and cx/, which have degree 1. We may associate with M a set ttm of half¬ 
edges of AM.n such that: (i) the midpoints Me and Mf have degree 1 , and (ii) every 
other vertex in AVn has even degree. Such a conhguration comprises a set of cycles 
together with a path between Me and Mf. The next lemma is immediate. 

Lemma 3.2. The two-edge correlation function of the 1-2 model satisfies 
(3-8) {eTe(Tf)n = ^ 

where the numerator Zn^e^^f is given in (2.7), and the parameters of the polygon 
model satisfy (3.7) and (3.5). 

3.4. The polygon model as a dimer model. We show next a one-to-one corre¬ 
spondence between polygon conhgurations on EI^ and dimer conhgurations on the 
corresponding Fisher graph of EI^. The Fisher graph is obtained from EI„ by 
replacing each vertex by a ‘Fisher triangle’ (comprising three ‘triangular edges’), as 
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illustrated in Figure 3.2. A dimer configuration (or perfect matching) is a set D of 
edges such that each vertex is incident to exactly one edge of D. 

Let TT be a polygon conhguration on ]HI„ (considered as a collection of edges). The 
local conhguration of vr at a black vertex n G is one of the four conhgurations at 
the top of Figure 3.2, and the corresponding local dimer conhguration is given in the 
lower line (a similar correspondence holds at white vertices). The construction may 
be expressed as follows. Each edge e of F„ is either triangular or is inherited from ]HI„ 
(that is, e is the central third of an edge of ]HI„). In the latter case, we place a dimer 
on e if and only if e ^ tt. Having applied this rule on the edges inherited from ]HI„, 
there is a unique allocation of dimers to the triangular edges that results in a dimer 
conhguration on F„. We write D = D{7r) for the resulting dimer conhguration, and 
note that the correspondence n ^ D is one-to-one. 



Figure 3.2. To each local polygon conhguration at a black vertex of 
]HI„, there corresponds a dimer conhguration on the Fisher graph F„. 
The situation at a white vertex is similar. In the leftmost conhguration, 
the local weight of the polygon conhguration is CbCc, and in the dimer 
conhguration A. 


By (2.2), the weight w{7i) is the product (over v G 14) of a local weight at v 
belonging to the set {eaCb, etCc, ecCa, 1}, where the particular value depends on the 
behavior of vr at n (see Figure 3.2 for an illustration of the four possibilities at a 
black vertex). We now assign weights to the edges of the Fisher graph F„ in such a 
way that the corresponding dimer conhguration has the same weight as vr. 

Each edge of a Fisher triangle has one of the types: vertical (denoted ‘V’), NE, or 
NW, according to its orientation. To each edge e of F„ lying in a Fisher triangle, we 
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allocate the weight: 

A if e is vertical, 

B if e is NE, 

C if e is NW, 

where A, B, C satisfy (3.6)-(3.7). The dimer partition function is given by 
(3.9) Z^{D) ;= 

D 

where D[s) C D is the set of dimers of type s. It is immediate, by inspection of 
Figure 3.2, that the correspondence tt y-)- D is weight-preserving, and hence 

= Z„(P). 

3.5. The spectral curve of the dimer model. We turn now to the spectral curve 
of the weighted dimer model on F„, for the background to which the reader is referred 
to [14]. The fundamental domain of is drawn in Figure 3.3, and the edges of F„ 
are oriented as in that hgure. It is easily checked that this orientation is ‘clockwise 
odd’, in the sense that any face of ]HI„, when traversed clockwise, contains an odd 
number of edges oriented in the corresponding direction. The fundamental domain 
has 6 vertices labelled 1, 2,..., 6, and its weighted adjacency matrix (or ‘Kasteleyn 
matrix’) is the 6x6 matrix W = {kij) with 

{ Wij if (f, j) is oriented from i to j, 

—Wij if (f, j) is oriented from j to f, 

0 if there is no edge between i and j, 

where the Wij are as indicated in Figure 3.3. From W we obtain a modified adjacency 
(or ‘modihed Kasteleyn’) matrix K as follows. 

We may consider the graph of Figure 3.3 as being embedded in a torus, that is, we 
identify the upper left boundary and the lower right boundary, and also the upper 
right boundary and the lower left boundary, as illustrated in the hgure by dashed 
lines. 

Let G C be non-zero. We orient each of the four boundaries of Figure 3.3 
(denoted by dashed lines) from their lower endpoint to their upper endpoint. The 
‘left’ and ‘right’ of an oriented portion of a boundary are as viewed by a person 
traversing in the given direction. 

Each edge {u, v) crossing a boundary corresponds to two entries in the weighted 
adjacency matrix, indexed {u,v) and {v,u). If the edge starting from u and ending 
at V crosses an upper-left/lower-right boundary from left to right (respectively, from 
right to left), we modify the adjacency matrix by multiplying the entry {u,v) by z 
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Figure 3.3. Weighted 1x1 fundamental domain of F„. The vertices 
are labelled 1 , 2 ,..., 6 , and the weights Wi^j and orientations are as 
indicated. The further weights are as indicated. 

(respectively, z~^). If the edge starting from u and ending at v crosses an upper- 
right/lower-left boundary from left to right (respectively, from right to left), in the 
modihed adjacency matrix, we multiply the entry by w (respectively, w~^). We 
modify the entry {v, u) in the same way. For a dehnitive interpretation of Figure 3.3, 
the reader is referred to the matrix following. 

The signs of these weights are chosen to reflect the orientations of the edges. The 
resulting modified adjacency matrix (or ‘modihed Kasteleyn matrix’) is 


/ -C B -I 0 0 \ 

C 0 -A 0 -z-^ 0 

-B A 0 0 0 -vj-^ 

1 0 0 0 -C B 

0 z 0 C 0 -A 

\ 0 0 w -B A 0 / 


The characteristic polynomial is given (using Mathematica or otherwise) by 


P{z, w) := det K 



(3.10) 
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By (3.6) and (2.8), 


P{z, ta) = 1 + + a^'y^ + /? 7 ^ + a'y{^ — 1) ( ta H- 


w 


+ 0/3(7^ - 1) + -^ + /37(a2 -!)[- + -]. 

\ z J \ z w/ 


The spectral curve is the zero locus of the characteristic polynomial, that is, the 
set of roots of P{z, w) = 0. It will be useful later to identify the intersection of the 
spectral curve with the unit torus = {{z,w) : lz\ = |ta| = 1 }. 

Let 


(3.11) 


U = a[3 + /dy + 70 — 1, 

V = -a/3 + /?7 + 70 ! + 1, 
S = a[3 — + 70 + 1, 

T = a[3 + /Sy — 70 + 1. 


Proposition 3.3. Let ea,eb,ec 7 ^ 0, so that a,(3,'y > 0. Either the spectral curve 
does not intersect the unit torus T^, or the intersection is a single real point of 
multiplicity 2. Moreover, the spectral curve intersects at a single real point if and 
only ifUVST = 0, where U,V,S,T are given by (3.11). 

Proof. The proof follows by a computation similar to those of the proofs of [11, Thm 
9] and [13, Lemma 3.2]. A number of details of the proof are very close to those of 
[11, 13] and are omitted. Instead, we highlight where differences arise. 

Let ea,eb,€c > 0. By (2.8) and (3.6)-(3.7), the map -0 : {A,B,C) 1 —>■ (0,0,7) 
is a bijection between (0, 00 )^ and itself. That P{z, w) > 0, for {z, w) G and 
(o, (3, 7 ) G (0, 00 )^, follows by the corresponding argument in the proofs of [11, Thm 
9] and [13, Lemma 3.2]. It holds in the same way that the intersection of P{z, w) = 0 
with can only be a single point of multiplicity 2 . 

We turn now to the four points when z,w = ±1. Note that 

P(l, 1) = (-1 + A^ + B^ + 

F(-l, -1) = {1-A^ + B^ + CY = S\ 

F(-l, 1) = {1 +A^ - B"^ + CY = T^, 

F(l, -1) = {1 + A^ + B^- CY = 


by (3.6). Since A,B,C 7 ^ 0, no more than one of the above four quantities can equal 
zero. □ 
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The condition UVST ^ 0 may be understood as follows. Let 7^ be given by (2.12), 
and note that 

(3.12) 72(a"\ /^) = l/7i(«, /5)- 

Proposition 3.4. Let Q!,/3,7 > 0 and let U,V,S,T satisfy (3.11). 

(a) We have that UV ST = 0 i/ and only i/7 G {71,72}- 

(b) The region i?sup of Theorem 2.5 is an open, connected subset of (0, cxd)^. 

(c) The region Rsnh is the disjoint union of four open, connected subsets of 
(0 , cxd)^, namely, 

^sub = {7 < 71} n {a/3 < 1}, i?2,,b = {7<7 i}n{a/ 3 >l}, 

RLh = {7 > 72} n {a < /3}, = {7 > 72} n {a > /3}. 

Proof. Part (a) follows by an elementary manipulation of (3.11). Part (b) holds since 
7i < 72 for all a, /3 > 0. Part (c) is a consequence of the facts that 71 = 0 when 
a + /3 = 0, and 72 = cxd when a — /3 = 0. □ 


4. Proof of Theorem 2.5 


4.1. Outline of proof. We begin with an outline of the main steps of the proof. 
In Section 4.2, the two-edge correlation of the polymer model on ]HI„ is expressed as 
the ratio of expressions involving Pfaffians of modihed adjacency matrices of dimer 
models on the graph of Section 3.4. The squares of such Pfaffian ratios are shown 
in Section 4.4 to converge to a certain determinant. This implies the existence of 
the limit M(e, /)^ of the two-edge correlation, and completes the proof of part (a) 
of the theorem. For parts (b) and (c), one considers the square M(e,/)^ as the 
determinant of a block Toeplitz matrix. By standard facts about Toeplitz matrices, 
the limit A := lim|e_/|_).oo Af(e,/)^ exists and is analytic except when the spectral 
curve intersects the unit torus. The remaining claims follow by Proposition 3.3. 


4.2. The order parameter in terms of PfafRans. By Remark 2.1, we shall as¬ 
sume without loss of generality that Ca, e?,, Cc > 0. Let £ be the path of AM.n connecting 
Me and Mf as in (2.11). To a conhguration vr G Llej we associate the conhgura- 
tion tt' := vr -|- £ G nP°h (with addition modulo 2). The correspondence vr -H- vr' is 
one-to-one between Llej and 11^°^. By considering the configurations contributing 
to Zn^e-i^f, we obtain by Lemma 3.2 that 


(4.1) 


MnieJ) 


Zn{P) 



Zn,e{P) 
Zn{P) ’ 


where Zn/{P) is the partition function of polygon conhgurations on AM.n with the 
weights of s-type half-edges along i changed from to e“^. 
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From the Fisher graph F„, we construct an augmented Fisher graph AWn by placing 
two further vertices on each non-triangular edge of F„, see Figure 4.1. We will con¬ 
struct a weight-preserving correspondence between polygon configurations on AM.^ 
and dimer configurations on AF„. 


✓ N 

X V 

✓ N 



Figure 4.1. The fundamental domain of AF„, which may be com¬ 
pared with Figure 3.3. 

We assign weights to the edges of AFn as follows. Each triangular edge of AFn is 
assigned weight 1. Each non-triangular s-type edge of the Fisher graph F„ is divided 
into three parts in AFn to which we refer as the left edge, the middle edge, and the 
right edge. The left edge and right edges are assigned weight while the middle 
edge is assigned weight 1. We shall identify the characteristic polynomial of this 
dimer model in the forthcoming Lemma 4.3. Let be the set of left and right non- 
triangular edges corresponding to half-edges in £, and let Vi be the set of vertices of 
edges in E^. 

There is a one-to-one correspondence between polygon conhgurations on AM.n and 
polygon conhgurations on EI^. The latter may be placed in one-to-one correspondence 
with dimer conhgurations on y4F„ as follows. Consider a polygon conhguration vr on 
]HI„. An edge e G is present in tt if and only if the corresponding middle edge 
of e is present in the corresponding dimer conhguration D = D{tt) on AF„. Once 
the states of middle edges of AF„ are determined, they generate a unique dimer 
conhguration on AFn- 

By consideration of the particular situations that can occur within a given funda¬ 
mental domain, one obtains that the correspondence is weight-preserving (up to a 
hxed factor), whence 
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where Zn{AD) is the partition function of the above dimer model on and Sg 

is the parameter corresponding to an edge with the type of g. A similar dimer 
interpretation is valid for Zn/{P), and thus we have 


(4.2) 


\ Z„AP) (TT -A Z'„(AD) 

z.(P) V Z„(P) I^II » ) ZMDY 


where Z'^{AD) is the partition function for dimer configurations on AF„, in which an 
edge of has weight eg (where g is the corresponding half-edge), and all the other 
left/right non-triangular edges have unchanged weights e~^. 

We assign a clockwise-odd orientation to the edges of AF„ as indicated in Figure 
4.1. The above dimer partition functions may be represented in terms of the Pfaffians 
of the weighted adjacency matrices corresponding to Zn{AD) and Z'^{AD). See 
[7, 8, 12, 21]. 

Recall that AF„ is a graph embedded in the n x n torus. Let and 'jy be two 
non-parallel homology generators of the torus, that is, and 'jy are cycles winding 
around the torus, neither of which may be obtained from the other by continuous 
movement on the torus. Moreover, we assume that y^, and y^ are paths in the dual 
graph that meet in a unique face and that cross disjoint edge-sets. For definiteness, 
we take 'jx (respectively, y^) to be the upper left (respectively, upper right) dashed 
cycles of the dual triangular lattice, as illustrated in Figure 4.2. We multiply the 
weights of all edges crossed by y^, (respectively, y^) by 2; or z~^ (respectively, w or 
w“^), according to their orientations. Note that i crosses neither y^. nor y^. 



Figure 4.2. Two cycles 'jx and y^ in the dual triangular graph of the 
toroidal graph ]HI„. The upper left and lower right sides of the diamond 
are identihed, and similarly for the other two sides. 
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Let Kn{z, w) be the weighted adjacency matrix of the original dimer model above, 
and let K'^{z,w) be that with the weights of s-type edges along i changed from 
to e^. 

If n is even, by (4.2) and results of [7, 12] and [18, Chap. IV], 

(4.3) 

Zn,e^f (n ,-A -Pf 1) + Pf 1) + Pf -1) + Pf "1) 

Zn{P) . j 2Z^{P) 

where 

(4.4) 2Z„(P) = -Pf iC„(l, 1) + Pf P„(-l, 1) + Pf P„(l, -1) + Pf Pn(-1, -1). 
The corresponding formula when n is odd is 

/tt _A Pf a';(i, 1) + Pf 1) + Pf -1) - Pfgg-i. -i) 

Z«(D (y ’ ) 2ZUP) 

as explained in the discussion of ‘crossing orientations’ of [19, pp. 2192-2193]. The 
ensuing argument is essentially identical in the two cases, and therefore we may 
assume without loss of generality that n is even. 


4.3. The limit as n —)• oo. In studying the limit of (4.3) as n —)■ cx), we shall require 
some facts about the asymptotic behaviour of the inverse matrix of iP„(6*, z/). We 
summarise these next. 

The graph may be regarded as n x n copies of the fundamental domain of 
Figure 4.1, with vertices labelled as in Figures 4.3-4.4. We index these by (p, q) with 
p,q = 1,2,... ,n, and let Dpq be the fundamental domain with index {p,q). Let 
P = {(P) ?) • D{p,q) O i ^ 0}, so that the cardinality of V depends only on \e — f\. 
Each Dp^q contains 12 vertices. For Vi,V 2 G Pi,i, we write Wi!-Dp2,g2 W2) 

for the {ui,U 2 ) entry of K~^, where Ui is the translate of Vi lying in Dp^^q.. 


Proposition 4.1. Let 6^, z/ G { — 1, !}■ We have that 

(4.5) lim K-^{e,iy){Dp^^q^,Vr; Dp,^^q^,Vs) 


1 

47J-2 




^P2 Plyjl2 


Vs,Vr 


dz dw 
iz iw 


where (pi, gi), (p2, ^2) G V, and r,s e {ai,bi : i = 1,2,3,4}, and 
denotes the {vs,Vr) entry of K^^{z,w). 


Proof. The limiting entries of K~^{6, z/) as n —)■ cx) can be computed explicitly using 
the arguments of [10, Thm 4.3] and [9, Sects 4.2-4.4], details of which are omitted 
here. □ 
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X 

✓ X 



Figure 4.3. The fundamental domain of with vertex-labels. 

\ 

\ 



N 

\ 


Figure 4.4. Part of the path I between two NW edges. 


Note that the right side of (4.5) does not depend on the values of 6^, z/ G {—1,1}. 


4.4. Representation of the PfafRan ratios. We return to the formulae (4.1) and 
(4.3)-(4.4) for the two-edge correlation. The matrices Kn{6,u) and K'^{9,v) are 
antisymmetric when E { — 1,1}. For 6,u E {—1,1}, 


(4.6) 



= det[Rr.K-\e,u) + l], 

where 


(4.7) 

Rn = K{e,u)-K^{e,iy). 


The following argument is similar to that of [12, Thm 4.2]. Let 
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and define the 4x4 block matrix 






s = a,b. 


Each half-edge of EI^ along i corresponds to an edge of namely, a left or right 
non-triangular edge. Moreover, the path i has a periodic structure in AM.^, each 
period of which consists of four edges of AEI^, namely, a NW half-edge, followed 
by two horizontal half-edges, followed by another NW half-edge. These four edges 
correspond to four non-triangular edges of A¥n with endpoints denoted Uai, 

L’aa, Uf,,, Vb^- See Figure 4.4. 

Let (p, q) G P. The 12 x 12 block of Rn with rows and columns labelled by the 
vertices in Dp q may be written as 


(Sa 0 0 \ 

(4.8) Rn{Dp^q, Dp^q) = j 0 —Sb 0 

\0 0 0 / 


Each entry in (4.8) is a 4 x 4 block, and the rows and columns are indexed by 
Va 2 , , Vai , L’ag, L’fea ) ) ^64 Wbg Wci, • • •, W • All other entries of Rn equal 0 . 

Owing to the special structure of Rn, the determinant of Sn '■= RnK~^{6, u) +1 is 
the same as that of a certain submatrix of Sn given as follows. From S'„, we retain all 
rows and columns indexed by translations (within V) of the and Vbj ■ Since each 
fundamental domain contains four such vertices of each type, the resulting submatrix 
Sn/ is square with dimension 8 11^1. By following the corresponding computations of 
[12, Sect. 4] and [18, Chap. VIII], we hnd that det S'„ = det 
Let Xg be the Vg x Vg block diagonal matrix with rows and columns indexed by 
vertices in Vg, and dehned as follows. Adopting a suitable ordering of Vg as above, 
the diagonal 2x2 blocks of Xg are Y (e* — e^^), where s depends on the type of the 
corresponding edge, and off-diagonal 2x2 blocks of Xg are 0. Note that 

(4.9) detX^ = n (^9 “ • 

Let i/)g be the submatrix of K~^{6, u) with rows and columns indexed by 

Vg. By Proposition 4.1, the limit 

(4.10) Kg := lim K-\e,iy)g 

n^oo 


exists and is independent of 6^, z/ G { — 1,1}. 

Proposition 4.2. The limit M{e,fY = lim„^oo Af„(e,/)^ exists and satisfies 


M{ejf 



(4.11) 
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Proof. Let 9,1^ E {—1,1}, and assume first that ea,eb 7^ 1. By (4.6)~(4.9) and the 
discussion before the proposition. 


detKf{9, u) 
det Kn{9, u) 


det[XiK-\9,u), + I] 

det[iL“^(6*, z/)^ + detX 

det[K-\e, iy)e + JJ 


On taking square roots, and noting that ^{9, z/)^ + ^ is antisymmetric. 


PiKjOx) 

PfKn{9,u) 


(-im Of+.c-‘i n • 

k ^gJ 


for some 7 that is independent of 9, v. 
By (4.3), 


2Z„,eo/ = (-l)''{-p(l,l)Pfi^n(l,l)+p(-l,l)PfX„(-l,l) 

+p(i, -i)pf x„(i, -1 )+p(-i, -i)pf x„(-i, -1)} j] [1 - ly 
where p(0, v) = Pf \K-^{9, v)i + Xj\ By (4.4) and (4.10), 

pt(A 7 ‘+jf,-‘) n ( 1 - 3 ) • 

g^Ei \ 5^. 

and (4.11) follows by (4.9) and (4.1). 

Assume next that Cq = e?, = 1. We have 

K'^{9,iy) = Kn{9,iy), for 0, z/= ±1. 

Since X7 = 0 in this case, we obtain (4.11) once again. If exactly one of e^, equals 
1, we obtain (4.11) as above. □ 


lim 

n^oo 


Zn{P) 


4.5. Proof of Theorem 2.5(a). As in [10, Thm 4.3] and [9, Sects 4.2-4.4], by 
Proposition 4.2, the limit M(e,/)^ = hm„_,.oo M„(e, /)^ exists and equals the deter¬ 
minant of a block Toeplitz matrix with dimension depending on |e — /|, and with 
symbol given by 

7(0 = ^ / 


( 4 . 12 ) 
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where T{(, 0) is the 8x8 matrix with rows and columns indexed by Vai,Va 2 ,Va 3 ,Va 4 ^, 
Vb^,Vb 2 ,Vb 3 ,Vb 4 (with rows and columns ordered differently) given by 






and Ag = 1 — e“^. See [23, 24, 25] and the references therein for accounts of Toeplitz 
matrices. 

One may write 


(4.13) 


[i^i ^{z,w)]ij 


Qij{z,w) 
P^{z, w) ’ 


where Qij{z, w) is a Laurent polynomial in z, w derived in terms of certain cofactors 
of Ki{z,w), and P"^{z,w) = det Ki{z,w) is the characteristic polynomial of the 
dimer model. 


Lemma 4.3. The characteristic polynomial P^ of the above dimer model on A¥n 
satisfies P^{z,w) = {eaebec)~'^P{z,w), where P{z,w) is the characteristic polynomial 
of (3.10). 

Proof. The characteristic polynomial P^ satisfies P^{z,w) = det Ki{z,w). Each 
term in the expansion of the determinant corresponds to an oriented loop configu¬ 
ration consisting of oriented cycles and doubled edges, with the property that each 
vertex has exactly two incident edges. It may be checked that there is a one-to-one 
correspondence between loop configurations on the two graphs of Figures 4.1 and 
3.3, by preserving the track of each cycle and adding doubled edges where necessary. 
The weights of a pair of corresponding loop configurations differ by a multiplicative 
factor of {ABCfi = {eaCbCcfi. □ 

By the above, the limit M{e,ffi exists whenever P^{z,w) has no zeros on the 
unit torus T^. By Lemma 4.3 and Proposition 3.3, the last occurs if and only if 
UVST 7^ 0. The proof of Theorem 2.5(a) is complete, and we turn towards parts 
(b) and (c). 


4.6. Proofs of Theorem 2.5(b, c). Consider an infinite block Toeplitz matrix J, 
viewed as the limit of an increasing sequence of finite truncated block Toeplitz ma¬ 
trices Jn- When the corresponding spectral curve does not intersect the unit torus, 
the existence of det J as the limit of det Jn is proved in [23, 24]. By Lemma 4.3 and 
Proposition 3.3, the spectral curve condition holds if and only if UV ST ^ 0. We 
deduce the existence of the limit 


(4.14) 


A(a,/3,7):=, lim 

\e-f\^oo 
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whenever UVST ^ 0. By Proposition 3.3, the function A is dehned on the domain 
D:= {0,ooY\{UVST = 0}. 

Lemma 4.4. Assume Q!,/3,7 > 0. The function A is an analytic function of the 
complex variables a, (3, 7 except when UV ST = 0, where U,V,S,T are given by 

(3.11). 

As noted after Theorem 2.5, A is singular when UV ST = 0. 

Proof. This holds as in the proofs of [12, Lemmas 4.4-4.7]. We consider A as the 
determinant of a block Toeplitz matrix, and use Widom’s formula (see [23, 24], and 
also [ 6 , Thm 8.7]) to evaluate this determinant. As in the proof of [ 6 , Thm 8.7], A 
can be non-analytic only if the spectral curve intersects the unit torus, which is to 
say (by Lemma 4.3 and Proposition 3.3) if UV ST = 0. □ 

The equation UV ST = 0 dehnes a surface in the hrst octant (0, 00 )^, whose 
complement is a union of hve open, connected components (see Proposition 3.4). 
By Lemma 4.4, A is analytic on each such component. It follows that, on any such 
component: either A = 0, or A is non-zero except possibly on a nowhere dense set. 
Let a,/9,7 > 0. By Proposition 3.4, UV ST 7 ^ 0 if and only if 

( 4 . 15 ) 7 e (0,7i) U (71,72) U (72,00), 

where the 7 * are given by ( 2 . 12 ). 

Proof of Theorem 2.5(b). By Proposition 3.4, UV ST 7 ^ 0 on the open, connected 
region i?sup- Therefore, A is analytic on i?sup- Hence, either A = 0 on i?sup, or A ^ 0 
on i?sup and the zero set Z := {r = (a,/?, 7 ) G i?sup : A(r) = 0 } is nowhere dense in 
i?sup- It therefore suffices to hnd (a,/?, 7 ) G i?sup such that A(q;,/ 3 , 7 ) 7 ^ 0. 

Consider the 1-2 model of Sections 3.1-3.3 with a = 6 > 0 and c > 4a. By (2.8), 
(3.5), and (3.6), the corresponding polygon model has parameters 

c — 2a U 

c + 2a’ ^ (c — 2 a)(c-|- 2 a) 

In this case, 72 = 00 and 7 G ( 71 , 72 )- 

By [ 6 , Thm 3.1], for almost every such c, the 1-2 model has non-zero long-range 
order. By Lemma 3.2, A(q;,/ 3 , 7 ) 7 ^ 0 for such c. 

Proof of Theorem 2.5(c). By Remark 2.2, when a ,/?,7 > 0 are sufficiently small, 
the two-edge correlation function M{e,f) of the polygon model equals the two-spin 
correlation function {o'eO'f) of a ferromagnetic Ising model on AH at high tempera¬ 
ture. Since the latter has zero long-range order, it follows that A = 0. Suppose, in 
addition, that a(3 < 1 and 7 < 71 . Since A is analytic on (in the notation of 
(3.13)), we deduce that A = 0 on Rl^^. We next extend this conclusion to R^^^ with 
k = 2,3,4. 
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Let (a,/ 3 , 7 ) G -Rgub- By (3.12), we have that a < 1 and 7 ^ < 71(0 ^,/5), so 
that (a“^,/ 3 , 7 “^) G .Rgub- By Theorem 2.3, 

A(a,/?, 7 ) = A(a“\/?, 7 “^) = 0. 

Therefore, A = 0 on 

Let (a,/ 9 , 7 ) G whence (a"^, 7 ) G i^g^b ^>7 (3.12). As above, 

A(a,/?, 7 ) = A(a"\/3"\7) = 0, 

whence A = 0 on i?s^b• The case of i?g^b can be deduced as was i?sub- 

Acknowledgements 

This work was supported in part by the Engineering and Physical Sciences Re¬ 
search Council under grant EP/I03372X/1. ZL acknowledges support from the Si¬ 
mons Foundation under grant 7)^:351813. The authors are grateful to two referees for 
their suggestions, which have improved the presentation of the work. 

References 

[ 1 ] R. J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic Press, London, 1982. 

[2] H. Duminil-Copin, R. Peled, W. Samotij, and Y. Spinka, Exponential decay of loop lengths in 
the loop 0(ji) model with large n, (2014), to appear. 

[3] M. E. Fisher, Statistical mechanics of dimers on a plane lattice, Phys. Rev. 124 (1961), 1664- 
1672. 

[4] G. R. Grimmett and S. Janson, Random even graphs. Electron. J. Gombin. 16 (2009), Paper 
R46, 19 pp. 

[5] G. R. Grimmett and Z. Li, The 1-2 model, (2015), http://arxiv.org/abs/1507.04109. 

[6] _, Critical surface of the 1-2 model, (2015), http://arxiv.org/abs/1506.08406. 

[7] P. W. Kasteleyn, The statistics of dimers on a lattice, I. The number of dimer arrangements 
on a guadratic lattice, Physica 27 (1961), 1209-1225. 

[8] _, Dimer statistics and phase transitions, J. Math. Phys. 4 (1963), 287-293. 

[9] R. Kenyon, Local statistics of lattice dimers, Ann. Inst. H. Poincare, Probab. Statist. 33 (1997), 
591-618. 

[10] R. Kenyon, A. Okounkov, and S. Sheffield, Dimers and amoebae, Ann. Math. 163 (2006), 
1019-1056. 

[11] Z. Li, Local statistics of realizable vertex models, Commun. Math. Phys. 304 (2011), 723-763. 

[12] _, Critical temperature of periodic Ising models, Gommun. Math. Phys. 315 (2012), 337- 

381. 

[13] _, 1-2 model, dimers and clusters, Electron. J. Probab. 19 (2014), 1-28. 

[14] _, Spectral curves of periodic Fisher graphs, J. Math. Phys. 55 (2014), Paper 123301, 25 

pp. 

[15] _, Uniqueness of the infinite homogeneous cluster in the 1-2 model. Electron. Gommun. 

Probab. 19 (2014), 1-8. 

[16] K. Y. Lin and F. Y. Wu, General vertex model on the honeycomb lattice: equivalence with an 
Ising model. Modern Phys. Lett. B 4 (1990), 311-316. 








CRITICAL SURFACE OF THE HEXAGONAL POLYGON MODEL 


25 


[17] M. Lis, Phase transition free regions in the Ising model via the Kac-Ward operator, Commun. 
Math. Phys. 331 (2014), 1071-1086. 

[18] B. McCoy and T. T. Wu, The Two-Dimensional Ising Model, Harvard University Press, Cam¬ 
bridge MA, 1973. 

[19] M. Schwartz and J. Bruck, Constrained codes as networks of relations, IEEE Trans. Inform. 
Th. 54 (2008), 2179-2195. 

[20] H. N. V. Temperley and M. E. Fisher, Dimer problem in statistical mechanics—an exact result, 
Philos. Mag. 6 (1961), 1061-1063. 

[21] G. Tester, Matchings in graphs on non-orientahle surfaces, J. Combin. Theory Ser. B 78 (2000), 
198-231. 

[22] B. L. van der Waerden, Die lange Reichweite der regelmdssigen Atomanordnung in Misch- 
kristallen, Zeit. Physik 118 (1941), 473-488. 

[23] H. Widom, On the limit of block Toeplitz determinants, Proc. Amer. Math. Soc. 50 (1975), 
167-173. 

[24] _, Asymptotic behavior of block Toeplitz matrices and determinants. II, Adv. Math. 21 

(1976), 1-29. 

[25] _, Toeplitz matrices and Toeplitz operators. Complex Analysis and its Applications, Lec¬ 

tures, Internat. Sem., Trieste (1975), vol. I, Internal. Atomic Energy Agency, Vienna, 1976, 
pp. 319-341. 

[26] X. N. Wu and F. Y. Wu, Exact results for lattice models with pair and triplet interactions, J. 
Phys. A 22 (1989), L1031-L1035. 

Statistical Laboratory, Centre for Mathematical Sciences, Cambridge Univer¬ 
sity, WiLBERFORCE ROAD, CAMBRIDGE CB3 OWB, UK 
E-mail address: g.r.grimmett@statslab.cam.ac.uk, 

URL: http://www.statslab.cam.ac.uk/~grg/ 

Department of Mathematics, University of Connecticut, Storrs, Connecticut 

06269-3009, USA 

E-mail address: zhongyang. liSuconn. edu 
URL: http://www.math.uconn.edu/~zhongyang/ 




